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THE  DISPERSION  OF  FIRE  SUPPRESSION  AGENTS  DISCHARGED  FROM  HIGH 
PRESSURE  VESSELS:  ESTABLISHING  INITIAL/BOUNDARY  CONDITIONS  FOR  THE 

FLOW  OUTSIDE  THE  VESSEL 


Leonard  Y.  Cooper 
Fire  Modeling  Group 
Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 


ABSTRACT 

This  work  reports  on  part  of  an  effort  to  study  the  dispersion  and  extinguishment  effectiveness  of 
Halon  and  Halon-altemative  fire  extinguishment  agents  discharged  from  N2-pressurized  vessels.  In 
the  systems  under  consideration,  as  the  agent  ©dts  from  the  vessel,  thermodynamic  and  fluid-dynamic 
instabilities  lead  to  flashing  and  break-up  of  the  agent  into  a  two-phase  droplet/gaseous  jet  mixture. 
This  occurs  in  a  transition  region  relatively  close  to  the  vessel  exit  oriflce/nozzle.  Downstream  of  this 
region  the  two-phase  agent  jet  then  mixes  with  the  ambient  air  environment  and  is  dispersed  in  the 
protected  space. 

A  mathematical  model  has  been  developed  previously  to  simulate  the  time-dependent  discharge  of 
the  agent  from  the  pressure  vessel  Using  the  output  of  this  model  and  thermodynamic  and  fluid- 
dynamic  considerations  of  the  phenomena  in  the  transition  section,  the  present  work  develops  a 
method  for  determining  a  set  of  initial/boundary  conditions  at  an  initial  section  of  the  jet,  downstream 
of  the  transition  region.  These  initial/boundaiy  conditions  are  in  a  form  that  can  be  used  to 
formulate  and  solve  the  problem  of  the  development  and  dispersal  of  the  ensuing  mixed  air/two- 
phase-agent  jet 

Example  applications  of  the  developed  methodology  are  presented.  These  are  for  agent  discharge 
from  a  half-liter  cylindrical  discharge  vessel  with  a  circular  discharge  nozzle/orifice  of  diameter 
0.019m.  Simulations  involve  discharge  of  the  vessel  when  it  is  half-filled  with  either  fteon  22  or 
Halon  1301  and  then  pressurized  with  N2  to  4137xlO^Pa  (600psi). 


Key  Words:  agents,  aircraft  fire  safety,  discharge,  fire  extinguishment,  fire  safety,  halon,  halon 
alternatives 


BACKGROUND  AND  PURPOSE 


The  Building  and  Fire  Research  Laboratory  (BFRL)  of  the  National  Institute  of  Standards  and 
(NIST^  is  carrying  out  a  Program  for  the  US  Air  Force  to  study  the  dispersion  and 
^inguishment  effectiveness  of  Halon  and  Halon-altemative  fire  extinguishment  agents  discharged 
from  N2-pressurized  vessels.  This  is  being  carried  out  in  support  of  an  effort  to  advance  the 
technology  of  fire  safety  in  US  Air  Force  aircraft. 

One  objective  of  the  program  is  to  measure  and/or  predict  the  dispersion  of  a  discharging  agent 
throughout  a  protected  space.  Tune-dependent  agent  concentrations  associated  with  a  particular 
agent  discharge  and  threat  scenario  would  be  used  as  a  basis  for  estimating  extinguishment 
effectiveness. 

This  work  focuses  on  one  aspect  of  the  problem  of  predicting  agent  dispersion  by  means  of 
mathematical/computer  modeling. 

The  strategy  to  predict  agent  dispersion  characteristics  involves  two  basic  mathematical  model 
components:  1)  a  model  to  simulate  the  time-dependent  discharge  of  the  agent  from  the  high- 
pressure  discharge  vessel;  and  2)  a  model  to  simulate  the  development  and  dispersal  of  the  ensuing 
mixed  au-/two-phase-agent  jet.  ® 

A  critical  element  in  exercising  any  component-2  model  is  to  establish  a  set  of  initiaVboundary 
condition  for  the  jet  This  must  be  at  a  location  of  the  jet  axis  near  the  exit  section  of  the  discharge 
vessel  ^e  imtiaVboundary  conditions  would  be  derived  from  the  predictions  of  the  component  1 
model.  The  pupxjse  of  this  work  is  to  develop  a  method  to  estimate  the  initial  conditions  from  the 
predictions  of  the  component-1  model  described  in  [1]. 


A  MODEL  OF  THE  DISCHARGE  PROCESS 

In  res^nse  to  the  above-mentioned  Program  objectives  a  component-l-type  mathematical  model  of 
agent  discharge  was  developed  m  [1].  The  model  is  designed  to  simulate  agent  discharge  for  test  or 
field-deployed  system  co^gurations  depicted  in  Figure  1.  This  involves  a  cylindrical  discharge  vessel 
with  a^ort,  circular,  exit  nozzle/orifice.  The  vessel  is  filled  with  pure  agent  (part  liquid  and  part 
gas).  TJe  agent  is  then  pressurized  with  N2  to  some  prescribed  value.  Liquid  agent  and  some 
dissolved  N2  K  m  the  lower  part  of  the  vessel.  Above  this  is  a  mixture  of  gaseous  agent  and  N,.  In 

the  ^  of  a  field-deploy^  system,  liquid-agent  discharge  is  initiated  by  actuation  of  an  explosivi  cap 
sealing  the  exit  nozzle/orifice.  ^  ^ 

to  simulate  field-deployed  system  discharges  under 
laboratory  conditions,  is  equipped  with  the  mdicated  high-pressure  N,  holding  tank.  This  is 
connected  to  the  discharge  vessel  via  an  orifice.  This  system  simulates  the  field-deployed  system  by 
using  a  rupture  diaphragm  over  the  exit  nozzle/orifice  (nominal  rupture  pressure  same  as  the  pre- 
charge  pr^ure  of  the  M-deployed  system)  rather  than  an  explosive  cap.  An  experimental  nm 
^gins  C  of  through-orifice  N2  flow  from  the  holding  tank.  Tie  vessel  is  pressurized  to 

the  point  of  diaphragm  rupture  and  this  is  immediately  foUowed  by  vessel  discharge  to  the 
atmc»Cnc  pressure  environment  Note  that  systems  of  interest  in  the  present  study  typically  involve 
mitial  discharge  vessel  pressures,  P^v,  of  several  tens  of  atmospheres,  where,  in  the  "upside  down" 


Figure-1  configuration  and  at  the  time  that  liquid  discharge  is  being  completed,  is  stQl  large 
enough  to  preclude  flashing  of  the  liquid  in  the  vessel. 

The  reader  is  referred  to  Reference  [1]  for  a  detailed  description  of  the  model  and  for  results  of 
simulations  carried  out  to  establish  an  experimental  design  and  procedure  which  would  be  expected 
to  closely  simulate  field-deployed  system  discharges. 

The  model  simulates  the  time-dependent  discharge  of  Figure-l-type  systems  up  to  the  time  that  the 
last  of  the  liquid  agent  is  ejected  from  the  vessel.  For  a  specified  set  of  geometric  system  parameters, 
agent  material,  and  initial  conditions,  key  results  of  a  model  simulation  are  the  rate  of  liquid  agent 
discharged  from  the  vessel  and  the  vessel  pressure.  This  work  develops  a  procedure  to  use  these 
results  along  with  other  considerations  to  estimate  initial/boundaiy  conditions  for  a  component-2-type 
model  simulation  of  the  development  and  dispersal  of  the  ensuing  mixed  air/two-phase-agent  jet 


MODELING  THE  EARLY  DEVELOPMENT  OF  THE  JET;  FROM  THE  EXIT  SECTION  OF  THE 
DISCHARGE  VESSEL  TO  A  NEARBY  SECTION  IN  A  STATE  OF  THERMODYNAMIC 
EQUILIBRIUM  AND  COMPLETED  DROPLET  FORMATION 

General  Considerations 

It  is  assumed  that  within  the  discharge  vessel  and  upstream  of  the  exit  nozzle/orifice  the  velocity  of 
the  liquid  agent  is  so  small  that  its  kinetic  energy  can  be  neglected  in  considerations  of  its 
therm^ynamic  state.  Therefore,  the  thermodynamic  state  upstream  of  the  nozzle/orifice  position  can 
be  estimated  by  the  pressure  in  the  discharge  vessel,  Pdv>  and  the  temperature  there  of  the  liquid 
agent,  Tj^y^.  Note  that  in  most  applications  it  is  racpected  that  Tpy^  is  well  approximated  by 
the  ambient  temperature  outside  the  vessel.  A  point  representing  the  state  of  the  liquid  agent 
in  the  vessel  is  indicated  in  the  sketched  pressure-enthalpy  (P,  h)  diagram  of  Figure  2. 

Because  of  the  relatively-short  nozzle\orifice  design  under  consideration,  the  time  interval  for  the 
liquid  to  approach  it,  pass  through  it,  and  enter  the  outside  environment  is  very  small.  For  conditions 
related  to  the  NIST  Program  this  is  of  the  order  of  lO'^s.  During  this  time  interval  the  pressure  of 
the  traversing  liquid  is  reduced  from  the  high  P^y  values  inside  the  vessel,  of  the  order  of  several 
tens  of  atmospheres,  to  pressures  of  the  order  of  P^mB’  atmosphere. 

As  discharging  liquid  enters  and  traverses  the  region  of  the  nozzle/orifice,  it  will  enter  and  move 
along  super-heated  metastable  thermodynamic  states.  While  dependent  on  its  changing  thermody¬ 
namic  state  in  the  vessel,  (P^^y,  Tj^y^y^,  and  on  its  particular  thermodynamic  properties,  as  the  agent 
material  penetrates  the  ambient  environment  it  approaches,  and  likely  achieves  the  P^mb  Pressure 
while  still  in  its  metastable  liquid  state.  Thus,  downstream  of  a  vena  contracta  the  liquid  agent 
develops  into  a  near-uniform-radius  liquid  jet,  ru,  which  can  be  described  following  traditional 
incompressible  fluid-dynamic  considerations.  See  Hgure  3.  For  the  relatively  large-vapor-pressure 
agent  materials  of  interest  here,  a  combination  of  fluid-dynamic  and  thermodynamic  instabilities  will 
then  lead  to  breakup  into  small  droplets  and  flashing  (Le.,  rapid  evaporation  to  a  two-phase  equilibri¬ 
um  thermodynamic  state)  of  the  metastable  liquid  jet  [2]. 

The  metastable  liquid  is  e}q)ected  to  move  through  thermodynamic  states  on  near-isentropic  paths. 
However,  as  discussed  in  [1],  for  a  given  agent  and  initial  state  conditions,  (P[)y,  T^y^y^,  it  is 
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possible  that  the  pressure  along  such  a  path  will  not  reach  P^mb  P^or  to  an  intersection  with  the 
agents  spinoidal  curve.  Moreover,  if  such  an  intersection  does  occur,  then  at  that  instant 
spontaneous  nucleation  in  the  liquid  is  to  be  expected  [3,  4].  It  is  conjectured  that  this  would  lead 
to  near-explosive  breakup  and  flashing  of  the  jet,  all  of  this  being  initiated  within  a  jet  penetration 
depth  (into  the  outside  environment)  of  the  order  of  a  single  nozzle/orifice  diameter. 

Whatever  the  themodynamic-state  path  of  the  discharging  metastable  liquid,  it  is  reasonable  to 
assume  that  a  relatively  short  distance  downstream  of  the  nozzle/orifice  a  thermodynamic  equilibrium 
state  and  an  end  to  the  droplet-formation  phenomenon  will  be  achieved.  It  is  also  reasonable  to 
assume  that,  as  the  agent  flows  further  downstream,  developing  as  a  mixed,  two-phase,  agent/air  jet, 
thennodynamic  equilibrium  will  be  maintained,  and  droplet  collision  and  agglomeration  will  not  play 
an  important  role  in  jet  dynamics. 


The  Initial  Section  of  the  Jet 

The  position  aloiig  the  jet  axis  where  toth  approximate  thermodynamic  equilibrium  and  completion 
of  droplet  formation  is  first  achieved  will  be  denoted  here  as  the  initial  position  of  the  jet  Associated 
wth  tlm  position  is  the  plane  section  normal  to  the  jet  axis,  called  the  initial  section  of  the  jet. 
Similarly,  the  thermodynamic  state  and  properties  of  the  jet  at  the  im'tial  position  will  be  denoted  as 
Its  initial  state  and  initial  properties.  The  region  of  the  flow  between  the  nozzle/orifice  discharge 
section  and  the  initial  section  will  be  referred  to  as  the  transition  region  of  the  jet. 

A  detailed  description  of  the  initial  state  and  properties  of  the  jet  would  be  extremely  complicated. 
Tnis  would  include  the  variation  across  the  initial  section  of  air/gaseous-agent  concentrations, 
temperature,  and  velocity,  and  of  the  size,  velocity,  and  temperature  distributions  of  the  liquid  agent 
droplets.  Note  also  that  in  the  present  application  all  of  these  variables  of  the  initial  jet  flow  field 
would  be  time-dependent  In  general,  an  accurate  detaUed  description  of  the  initial  state  of  the  jet, 
whether  determined  by  theoretical  or  experimental  means,  is  either  impractical  or  beyond  the  current 
state  of  technology. 


Although  it  is  not  possible  to  provide  a  detailed  description  of  the  initial  state  of  the  jet,  it  is 
reasonable  to  aq)ect  that  an  achievable  approximate  description,  suitable  for  use  in  a  component-2- 
type  model,  would  yield  good  approximations  of  the  dispersing-agent-flow-field  problem.  It  is  the 
purpose  of  this  work  to  provide  such  an  approximate  description. 


^sumptions  of  the  Phenomena  in  the  Tkansitlon  Region:  Fkom  the  Nozzle/Orifice  to  the  Initial 
Section 

The  djmamic  processes  that  uutiate  breakup  and  flashing  of  the  metastable  Uquid  jet  and  bring  it  to 
Its  imtial  state  are  not  completely  understood.  At  one  extreme,  it  is  possible  that  the  above- 
mentioned  instabffities  lead  to  violent  breakup  and  flashing  of  the  hquid  jet  immediately  upon  leaving 
the  exit  nozzle/orifi^  Such  behavior  was  observed,  for  example,  in  some  experiments  of  151.  Also 
toward  the  completion  of  the  liquid  discharge,  this  was  also  observed  in  recent  experiments  at  NIST 
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The  metastable  state  can  also  persist  even  as  the  liquid  jet  penetrates  relatively  deeply  into  the 
ambient  environment.  For  example,  Reference  [5]  describes  a  liquid  water  jet  of  diameter  D  = 
7.6(10‘^)m,  initially  at  9.03xl0^Pa  (131psi)  and  415K,  which  penetrated  up  to  0.0254m,  corresponding 
to  L/D  »  32,  before  initiation  of  violent  breakup/flashing.  Observed  in  the  recent  NIST  experiments 
[6]  was  a  liquid  Freon  22  jet,  D  «  0.02m,  initially  at  41.37xlO^Pa  (600psi)  and  294K,  which  penetrated 
a  distance  of  approximately  0.3m,  corresponding  to  L/D  «  15,  before  initiation  of  violent  break¬ 
up/flashing.  In  this  regard  it  is  of  interest  to  note  that  in  the  literature  of  the  discharge  of  flashing 
liquids  from  high  pressure  vessels  there  is  a  degree  of  conventional  wisdom  that  analyses  of  the  two- 
phase  flow  phenomena,  based  on  thermodynamic  equilibrium  states,  tend  to  yield  good  results  when 
the  nozzle/piping  from  the  inside  to  the  outside  of  a  high  pressure  discharge  vessel  is  at  least  0.1m 


The  initial  section  will  be  some  incremental  distance  downstream  of  the  penetration  distance  of  the 
metastable  liquid  jet  The  studies  of  [5]  and  the  recent  experiments  at  NIST  [6]  suggest  that  this 
incremental  distance  will  typically  not  exc^  the  order  of  the  penetration  distance  itself.  Whatever 
the  distance,  it  is  assumed  here  that  the  total  length  of  the  transition  region  is  small  enough, 
compared  to  the  characteristic  length  of  the  overall  agent  dispersal  problem  of  interest,  that  the 
processes  within  the  transition  region  can  be  treated  as  quasi-steady. 

fa  the  transition  region  it  is  assumed  that  momentum  transfer  and  any  net  work  or  heat  transfer 
interactions  between  the  penetrating  agent  jet  and  the  ambient  air  environment  are  negligible.  Thus, 
consistent  with  the  analji^is  of  [8]  for  steady,  two-phase,  air/evaporating-liquid  jets,  it  is  assumed  that 
the  jet  at  the  initial  section  can  be  approximated  as  consisting  only  of  a  two-phase  mixture  of  agent, 
with  no  entrained  air.  Finally,  in  the  transition  section  it  is  assumed  that  increases  in  the  energy  of 
the  liquid  due  to  mechanical  work  by  surface  tension  forces  to  form  the  initial-section  droplets  are 
negligible  compared  to  the  energy  transfers  involved  in  phase-change  processes. 


The  Initial  State  -  A  Summary  of  Assumptions  and  Approximations  and  Their  Implications 

Regarding  ^  estimate  of  the  imtial  state,  the  assumptions  of  the  previous  paragraph,  some  additional 
approximations  and  assumptions,  and  their  implications  are  now  summarized. 


At.  the  initial  section.  The  jet  is  defined  by  a  characteristic  initial  radius  ty.  The  region  r  S  n 
consists  of  a  pure  mixture  of  saturated-agent  vapor  and  liquid  (droplets)  in  thermodynamic 
equOibrium  at  temperature,  Tj,  and  pressure,  Pj,  where 


“  ^AMB»  (1) 

“  ^SAT  “  r  s  ty  (2) 


and  where 
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^SAT  ~  ^  ^  ^CR 


(3) 


is  the  saturation  temperature  of  the  agent,  a  specified  function  of  pressure  P,  where  P  can  not  exceed 
the  critical  pressure,  P The  liquid-agent  component  of  the  jet  consists  of  droplets  of  uniform- 
radius  rj  u  and  these  are  uniformly  dispersed  over  the  region  r  £  rj.  The  agent  droplets  and  gas  have 
the  same  uniform  axial  jet  velocity,  Vj^. 


^  (4) 


where  Vj  =  V,(r)  is  the  velocity  distribution  of  the  jet  at  the  initial  section.  There  is  no  agent  in  r 

>  rj. 

The  region  r  >  rj  consists  of  a  uniform  ambient-temperature/pressure  environment  fPA»*n  T»»,t, 
of  quiescent  air. 


^AMb5  Tj  —  —  0;  r  >  rj 

There  is  no  air  in  r  S  rj. 


(5) 


Conservation  of  ma?8s; 


Pl^  =  l/[x,/yO^AG  +  (1  -  Xl)/Pl>d 


where pj*,  />j^  and pj^Q  are  the  average  density  of  the  agent,  the  density  of  the  liquid  agent 
(droplets),  Md  the  density  of  the  agent  gas  at  the  initial  section,  respectively;  x,  is  the  quality  of  the 
^o-p  flow  at  the  initial  section,  Le.,  the  fraction  of  a  mass  of  agent  which  is  saturated  gas; 

^  ,  instantaneous  mass  of  liquid  agent  in  the  discharge  vessel;  and  the  right  side  of  Eq' 
(6),  detenmned  previously  from,  e.g.,  the  discharge  model  of  [1],  is  assumed  to  be  specified. 


Conservation  of  energy  fand  using  Fk|. 


•»DV,AL  ~  *»I,AL  +  *lhALG(TsAr)  + 

^algCTsat)  =  -  hj^ 
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where  hp^AL  speciGc  enthalpy  of  the  agent  liquid  in  the  discharge  vessel  and  at 

the  initial  section,  respectively;  hj^o  is  the  specific  enthalpy  of  the  agent  gas  at  the  initial  section; 
and  h^Q  is  the  heat  of  vaporization  of  the  agent  at  the  temperature  TgAp  In  Eq,  (9), 


^ALG  ~  ^  Tcr  (10) 

is  a  specified  function  of  temperature,  X  where  T  can  not  exceed  the  critical  temperature, 

A  point  representing  the  state  of  the  agent  at  the  initial  section  is  presented  in  the  (P,  h)  diagram 
of  Figure  2.  Consistant  with  Eq.  (8),  this  reflects  the  fact  that  some  of  the  original  enthalpy  of  the 
material  in  its  high  pressure  "rest"  state  was  ecchanged  for  its  now-significant  kinetic  energy. 

The  sketch  of  the  transition  region  and  of  the  initial  properties  as  presented  in  Figure  3  is  consistent 
with  the  all  the  above  assumptions/approximations. 


Droplet  size.  It  is  assumed  that  rjo  is  specified.  Based  on  [5]  it  is  expected  that  this  will  be  in  the 
range  lOpm  -  100pm.  Photograpfiic  data  during  future  NIST  experiments  should  provide  sharper 
estimates  of  this  parameter.  Also,  the  sensitivity  of  the  downstream  agent  concentrations  to  variations 
of  this  parameter  should  be  evaluated  with  component-2  model  simulations. 


Initial  jet  radius.  It  is  assumed  that  ty  is  specified.  Based  on  photographic  data  in  [5]  and  on 
photographic  data  alreacty  acquired  during  NIST  experiments,  rj/ru  seems  to  be  of  the  order  of  10  - 
20.  Thus,  in  (what  is  reasonable  to  construe  visually  in  high  speed  photographs  as)  a  region  of 
breakup  and  flashing  of  the  liquid  jet  to  an  initial  equilibrium  state,  the  radius  of  the  jet  increases  by 
a  factor  of  10  -  20.  Again,  the  sensitivity  of  the  downstream  agent  concentrations  to  variations  of  this 
ratio  should  be  evaluated  with  component-2  model  simulations.  One  would  hope  that  such  sensitivity 
is  not  great 


^rameters  defining  the  Initial  state.  Consistent  with  the  above,  the  specified  parameters  of  the 
initial  state  are  Pj,  according  to  Eq.  (1),  r^j,  and  TjA'lj.  Also,  for  a  particular  agent  of  interest,  Tj 
is  specified  by  (2)  and  (3).  Required  to  complete  the  description  of  the  initial  state  are  Pu^  ^LA» 


Completing  the  Description  of  the  Initial  State 

MQdelin£  the  lionid  agent  as  incompressible.  For  problems  of  interest  here,  agent  liquid 
temperatures  are  expected  to  be  small  enough  compared  to  T^r  and  vary  in  a  small  enough  range 
to  permit  accurate  modeling  of  the  liquid  as  being  incompressible,  with  density  pp^  and  as  having 
a  constant  specific  heat,  C^.  Thus, 
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~  Pal 


(11) 


The  values  of  pj^  and  Cal  will  be  taken  as 


Pal  /^alC’I’al);  C^  -  C^(T^;  +  T^y^/Z  (12) 


where  is  taken  to  be  the  average  temperature  of  the  liquid  in  the  discharge  vessel  and  at  the 
initial  state,  and  where and  C^(T)  are  specified  functions  of  T 

From  the  assumption  of  incompressibility  it  follows  that 

l^DYAL  *  1*I,AL  “  C^(Tj)y^  -  TjyQ-)  +  (Pdv  •  PaMb)4'’aL  (13) 


Modeline  the  agent  gas  as  a  perfect  gas.  The  agent  gas  is  modeled  as  a  perfect  gas. 


PlA^  1’i/(1^AgTi)  ~  ^ambI^aq^sa^ 


(14) 


where  is  the  gas  constant 


T^e  solution  for  Vj^^^and^.  As  discussed  above,  it  is  expected  that  in  most  instances  the  liquid  jet 
will  achieve  P^mb  I*®  breakup/flashing.  Then,  Bernoulli’s  equation  leads  to 

dMoYAL/dt  =  ‘>^uW2(PdV  -  PAMB)//>Alj^  (15) 


where  Tjj  is  determined  from  the  radius  of  the  exit  nozzle/orifice 
coefficient,  Cp. 


,  rj^,  and  its  oriGce  discharge 


(^u^n)^  -  Cj} 

Substituting  Eqs.  (9)  and  (13)  into  Eq.  (8)  lead  to 

Vj^  a  =  Cal(Tdv^  -  Tg^)  +  (Pjjy .  Pamb)//>al  -  xihALG(TsAr) 


(16) 


(17) 
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Substituting  Eqs,  (7)  and  (15)  into  Eq.  (6)  leads  to 


~  (*‘u/*'i)'*[(Pdv  -  PambV/®alI[(V®)(1  -  «)]^  (18) 

®  =  PiaqIPa^  (19) 

Comparing  Eqs.  (18)  and  (17)  leads  to 

+  A2X  -  (Aj  +  1)  =  0  (20) 

where 


X  =  (Xi/e)(l  -  £)(ru/r,)2;  Xj  =  Xe/[(1  -  6)(rLj/ri)2] 

(21) 

h  -  Pal^ajJJhwj^  -  TsatVCPov  *  ^amb) 

(22) 

A2  =  [e/(l  -  OIPal^algCTsatMC^dv  ‘  ^amb)(’'u^i)^1 

(23) 

The  solution  to  Eq.  (20)  is 


X  =  -  (A2/2)  +  +  Aj  +  ifP 


(24) 


With  (21)  and  other  previous  definitions,  Eq.  (24)  provides  the  solution  for  x,.  This  would  be 

used,  in  turn,  to  obtain  Vj^  from  Eqs.  (17)  or  (18),  thereby  completing  the  description  of  the  initial 
state. 


^i,A  ~  2X[(Pj3y  -  (17’) 

A  sketch  of  the  velocity  distribution  and  the  mass  flux  distributions  of  gaseous  and  liquid  agent  at  the 
initial  section  are  presented  in  Figure  4. 


ALGORITHM  FOR  DETERMINING  THE  INITIAL  STATE 
The  above  solution  can  be  implemented  with  the  following  algorithm- 
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1.  Specify  agent  and  its  thermodynamic  properties,  including  Tg^^P)  and  h^G(T),  according 
to  Eqs.  (3)  and  (10),  respectively,  and/?^(T^,  Cal(T). 

2.  Specify  rj/rjj  (expected  range  10  -  20). 

3.  Specify  and  T^v^.  (In  most  cases  expect  to  specify  =  T^^.) 

4.  Specify  Tj  =  T^mb  for  r  >  fi*  Compute  from  Eq.  (3)  and  then  specify  Tj  = 

"^SAT  for  0  ^  r  S  rj. 

5.  Compute  hy^G(Ts^)  from  Eq.  (10). 

6.  Compute  T^l  according  to  Eq.  ( 12)  and  then  and  C^^  according  to  Eq.  (12)  and  PalCI) 
and  Cal(T),  respectively. 

7.  Compute  according  to  Eq.  (14)  and  then  e  according  to  Eq.  (19). 

8.  Use  a  component-1  model,  e.g.,  that  of  [1],  to  find  the  time-dependent  values  of  P^y  for  a 
particular  Figure-l-configuration  of  interest 

9.  For  a  particular  value  of  P^y  compute  Aj  and  X2  from  Eqs,  (22)  and  (23),  and  then  X  from 
Eq.  (24). 

10.  Fmd  xj  from  Eq.  (21)  and  Vj^  from  Eq.  (17’), 


EXAMPLE  APPUCATION  OF  THE  ALGORITHM 

Conditions  for  Example  Discharges  and  Temperature-Dependent  Agent  Property  Functions 

The  algorithm  was  applied  in  four  example  discharge  simulations.  These  involve  the  agents  R-eon 
22  (CHCIF2)  and  Halon  1301  (CBrFg)  and  two  different  T]  values  corresponding  to  rj  =  IOtjj  and 
20rLj.  The  examples  consider  time-dependent  agent  discharges  from  a  Figure-1  configuration  with 
no  holding  tank.  The  geometric  parameters,  orifice  coefficient,  and  initial  conditions  assumed  in  the 
simulated  example  discharge  calculations  are  presented  in  Table  1.  These  are  the  same  parameters 
used  in  the  example  calculations  of  [1].  During  the  discharge,  the  time-dependent  value  of  Pr,y  was 
computed  using  the  component-l-type  model  of  [1]. 

The  material  property  functions  for  Freon  22  (CHCIF2)  and  Halon  1301  (CBrFg)  which  were  used 
in  the  example  calculations  are  from  [9].  These  are  of  the  form 


Molecular  weight,  M:  M  kg  of  agent  =  1  kg-mole  of  agent  (25) 

P(TsAr)/Pa  =  exp[A  +  B/(Ts^)  +  Cln(TsA3yR)  +  DCrsAjVR)^];  Tj  s  Tsat  ^  ^2 

(26) 
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Complete  descriptions  of  these  functions  for  the  two  agents  are  presented  in  the  APPENDIX. 

Tsat  obtained  as  the  root  of  Eq.  (26)  where,  according  to  Eq.  (2)  or  item  4  of  the  algorithm,  the 
left  hand  side  was  specified  as  Py^.  Then  T^l  was  calculated  using  Eq.  (12).  Finally,  hALG(TsAr). 
Pal(JaD^  and  C^(Tal)  were  calculated  fi-om  Eqs.  (27),  (28),  and  (29),  respectively.  The  results 
are  presented  in  Table  2. 


Results  of  the  Example  Calculations 

Die  results  of  the  example  calculations  are  tabulated  in  Dble  3  (R-eon  22,  r.  =  lOri,  =  lOCn^^rw 
=  0.07378m),  Dble  4  (Halon  1301,  r,  =  IOtu  =  lOCp^j^  =  0.07378m),  Tkble  5  (Freon  22,  r,  = 
20ru  =  20Ci5^j^  =  0.1476m),  and  Tkble  6  (Halon  1301,  Tj  =  20rLj  =  20CD^rj,  =  0.1476m. 
Included  in  the  tables  are  computed  values  of  the  variable  5,  which  is  the  dimensionless  mass  of  liquid 
agent,  or  dimensionless  elevation  of  the  liquid  agent  in  the  discharge  vessel,  Z. 


^  -  Mdv,al^dv,al(*  -  0)  *  Z/Z(t  -  0)  (30) 


The  tabulated  results  of  Tkble  3  for  the  discharge  of  Freon  22  (tj  =  lOr^  =  0.07378m)  are  plotted 
in  Figures  5  -  8.  These  include  plots  of  the  time-dependent  values  of  Ppy  (Figure  5),  Vj.  (Figure 
6),  5  and  Xj  (Figure  7),  and  dMj^/dt  and  dMj^o/dt  (Figure  8). 

A  stucfy  of  the  results  leads  to  the  following  observations: 

1.  For  both  Freon  22  and  Halon  1301,  the  times  to  complete  the  discharge  of  liquid  agent  from 
the  pressure  vessel  are  almost  identical  at  0.024  ±  0.001s  (Tkbles  3  -  6). 

Z  At  the  initial  section,  the  quality  of  the  two-phase  agent  flow,  Xj,  is  significantly  different  for 
the  two  agents  considered,  but  for  a  given  agent  Xj  is  relativefy  independent  of  t  and  rj,  for 
the  two  initial-section  radii  values  considered  (rj  =  0.07378m  and  0.1476m).  Thus,  for  the 
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cases  considered,  x,  =  0.31  ±  0.01  for  Freon  22  (Tables  3  and  5)  and  x,  =  0.54  ±  002  for 
Halon  1301  (Tables  4  and  6).  ^ 

At  the  initial  section,  for  a  given  agent,  and  for  the  two  initial-section  radii  values  considered, 
the  mass  flow  rates  of  liquid  or  gaseous  agent  are  almost  identical  at  any  particular  time 
during  the  discharge.  [R»r  example:  for  Freon  22,  approximately  half-way  through  the 
=  ®*°1252s,  values  of  dMj^/dt  for  rj  =  0.07378m  and  rj  =  0.1476m  are 
8.59kg/s  (Table  3)  and  8.49kg/s  (Table  5),  respectiveljr,  and  for  Halon  1301,  approximately 
half-way  through  the  discharge  at  t  =  0.01222s,  values  for  dMj^/dt  for  rj  =  0.07378m  and 
ri  =  0.1476m  are  6.95kg/s  (Table  4)  and  8.59kg/s  (Table  6),  respectively.]  However,  over  the 
Murse  of  the  discharge,  these  mass  flow  rates  are  reduced  to  approximately  two  thirds  of  their 
mitial  value.  [For  example:  for  Freon  22  and  rj  =  0.07378,  the  values  of  dM,  Wdt  at  t  =  0 
and  t  =  02326s  (near  the  end  of  the  liquid  discharge)  are  11.78kg/s  and  7.34kg;S;>espectively 
^able  3  and  Figure  5);  and  for  Halon  1301,  and  r,  =  0.07378,  the  values  of  dM,^  /dt  at  t 
“  0  and  t  -  0.2443s  (near  the  end  of  the  liquid  discharge)  are  9.26kg/s  and  5  88kg/s 
respective^  (Table  4).]  ^  ’ 
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NOMENCLATURE 


A,B,C,D  constants,  Eqs.  (A2)  -  (A5) 

constant  specific  heat  of  liquid  agent 
orifice  coefficient 
diameter  of  a  liquid  jet 
Eqs.  (28)  and  (A4) 
specific  enthalpy 

specific  heat  of  vaporization  of  the  agent 
h  of  agent  liquid  in  the  discharge  vessel 
h  of  agent  gas  (liquid)  at  the  initial  section 
length  of  liquid  jet 
molecular  weight 

(dMj^/dt)  mass  flow  rate  of  agent  gas  (liquid)  at  the  initial  section 


^AL 

Cd 

D 

F 

h 

*^ALG 

^DV,AL 

L 

M 


M, 


DV,AL 


■AMB 


CR 


'■AG 


^03 


mass  of  liquid  agent  in  the  discharge  vessel 
pressure 

p  of  the  ambient  environment 

critical  pressure 

p  in  the  discharge  vessel 

p  at  the  initial  section 

gas  constant  for  the  agent  gas 

radius 

r  of  two-phase  agent  jet  at  the  initial  section 
r  of  agent  droplets  at  initial  section 
r  of  the  liquid  jet  outside  the  discharge  vessel 


14 


T 

Tal 

Tamb 

"^DYAL 

Ti 

^SAT 

Ti,T2 


t, 


BURST 


z 


z 


Zl 

e 

Ai,  A2 
PlA 

PlMi  ^IaD 

Pal 

C 


r  of  the  discharge  vessel’s  exit  nozzle/orifice 

absolute  temperature 

T  of  liquid  agent 

T  of  ambient  environment 

T  of  liquid  agent  in  the  discharge  vessel 

T  at  the  initial  section 

saturation  temperature  of  the  agent 

constant  temperatures,  Eqs.  (A2)  -  (A5) 

time 

t  when  the  diaphragm  bursts 
velocity  at  the  initial  section 
Vj  of  agent 
Vj  of  air 
Eq.  (21) 

quality  of  the  two-phase  agent  flow  at  the  initial  section 

distance  of  the  liquid/gas  agent  interface  above  the  exit  nozzle/orifice 

distance  along  the  jet  axis  from  the  exit  orifice/nozzle 

z  at  the  initial  section 

Eq.  (19) 

Eqs.  (21)  and  (22) 

average  density  of  the  agent  jet  at  the  inital  section 
density  of  the  agent  liquid  (gas)  at  the  initial  section 
density  of  liquid  agent 

dimensionless  elevation  of  the  liquid/gas  interface,  Eq.  (30) 
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APPENDIX:  MATERIAL  PROPERTY  FUNCTIONS  FOR  FREON  22  (CHCIFj)  AND  HALON 

1301  (CBrFj) 


The  following  material  property  functions  for  Freon  22  (CHCIF2)  and  Halon  1301  (CBrFj)  were  used 
in  the  example  calculations  [9]: 


Molecular  weight  =  M  kg  of  agent  =  kg-mole  of  agent 


(Al) 


M 


Freon  22  ('CHC1F2'>  Halon  1301  fCErF^'t 

86.468  148.910 


P(TsAr)/Pa  =  exp[A  +  B/(Ts^)  +  Cln(Ts^)  +  DCTsAiyR)^]:  T^  <  Tg^  ^  T2 

(A2) 


A 

B 

C 

D 

Ti 

T2 


Fi-eon  22  fCHClF^l 

9.7882(10) 

-  4.7735(10-^) 

-  1.2317(10) 
2.4173(10'^) 
115.73K 
369.30K 


Halon  1301  rCBrFjl 

6.2780(10) 

-  3.3770(10^) 

-  6.7461 
1.4407(10-^) 

105.15K 

340.15K 
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hALG(TsArV(J/kg-mole) 


“  ^ALG(TsArV(J'^g) 

=  A(1  -  Tsy^j/To^)®;  Tj  £  Tg^  <  T2 


(A3) 


A 

B 


/>AL(T)/(kg-mole/m^) 


A 

B 

C 

D 

Ti  = 

T2  = 


Freon  22  rCHaP.'l 

2.9617(10^) 

3.8600(10-1) 

369.30K 

115.73K 

369.30K 


Reon  22  rCHClF^'l 

1.5956 

Z6605(10*l) 

3.6930(102) 

Z8123(10*l) 

115.73K 

369.30K 


Halon  1301  rCBrF^,') 

2.5070(10’^) 

3.5700(10-1) 

340.15K 

105.15K 

340.15K 


Halon  1301  rCBrF^^ 

13435 

2.6995(10-1) 

3.4015(102) 

18006(10-1) 

105.15K 

340.15K 


=  /OAL(T)/(Mkg/m2) 

=  A/B^;  F  =  {1  +  [1  -  (T/R)/ql^};  Tj  <  T  <  T2 


(A4) 
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CAL(W/(kg-mole-K)]  =  MC^LCW/Ckg-K)]  (A5) 

=  A  +  B(T/R)  +  C(T/R)2  +  D(T/R)^  <  T  <  T2 


Freon  22  ('CHCIF.'l 

Halon  1301  fCBrFA 

A 

9.4660(10^) 

1.1421(10^) 

B 

3.9300(10^) 

-  1.3150(103) 

C 

-  8.1740(10'^) 

6.1200(10-1) 

D 

2.7910(10-3) 

0 

Ti  = 

115.73K 

193.15K 

T2  = 

369.30K 

298.15K 
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Ibble  1. 


Geometric  parameters,  orifice  coefficient,  and  initial  conditions  for  example 
discharges  from  a  Figure-l-type  configuration  with  no  holding  tank. 


Length  of  discharge  vessel  =  Zjjy  = 

0.2546m 

Cross-section  area  of  discharge  vessel  =  Apy  = 

0.1963xl0-2m2 

Volume  of  discharge  vessel  =  V^y  = 

0.50Qxl0'^m^ 

Area  of  exit  nozzle/orifice  =  Aq  = 

0.285Qxl0-V 

Radius  of  exit  nozzle/orifice  =  Tq  = 

0.9525xl0'^m 

Discharge  coefficient  of  exit  nozzle/orifice  =  Cjj  = 

0.6 

Initial  volume  of  agent  liquid  in  discharge  vessel  =  Mjjy^  j  = 

0.250xl0-^m^ 

Initial  temperature  of  gas  in  the  discharge  vessel  =  Tjjyj  = 

294K 

Initial  pressure  in  discharge  vessel  =  P^yj  = 

41.36xlO^Pa  (600psi) 

Temperature  of  agent  liquid  in  discharge  vessel  =  Tpy^  = 

294K 

Table  2.  Constant  properties  of  the  initial  state  for  the  example  calculations. 


Freon  22  rCHClF,'> 

Halon  1301  rCBrF^'i 

^SAT 

232K 

215K 

Tal 

263K 

254.5K 

^algCTsat) 

0.234xl0<^J/kg 

0.118xl0*^J/kg 

/^al(Tal) 

0.132xl(f‘kg/m^ 

0.1g0xl0^kg/m^ 

^alCTal) 

0.115xl0^J/(kg-K) 

0.809xl0^J/(kg-K) 
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I^ble  3. 


Initial  state  of  a  Freon  22  (CHCIF2)  jet  from  Figure-l-type  vessel  with  characteris¬ 
tics  of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  r,  =  lOrr  t  = 
=  0.07378m.  ^ 


Pal  =  1320kg/m3 
PViG  =  3.58kg/m^ 
T^a  =  232K 


t 

[S] 

^DV 

[Pa] 

c 

Vu 

[m/s] 

*1 

[kg/s]  [kg/s] 

0 

0.4135E+07 

1.0000 

122.73 

03022 

5.10 

11.78 

0.00179 

0J663E+07 

0.9031 

115.46 

03025 

4.80 

11.07 

0.00358 

0J299E+07 

0.8117 

109.49 

03028 

4.55 

10.48 

0.00537 

0.3008E+07 

0.7248 

104.46 

03030 

434 

9.99 

0.00716 

0.2769E+07 

0.6417 

100.14 

03032 

4.16 

9.57 

0.00895 

0.2570E+07 

0.5620 

9637 

03033 

4.01 

9.20 

0.01073 

0.2401E+07 

0.4851 

93.05 

03034 

3.87 

8.88 

0.01252 

0.2255E+07 

0.4108 

90.08 

03035 

3.75 

8.59 

0.01431 

0.2128E+07 

03389 

87.41 

03036 

3.63 

833 

0.01610 

0.2016E-h07 

03690 

84.99 

03037 

3.53 

8.10 

0.01789 

0.1917E-h07 

03010 

82.78 

03038 

3.44 

7.89 

0.01968 

0.1829E+07 

0.1347 

80.76 

03039 

336 

7.69 

0.02147 

0.1749E+07 

0.0701 

78.88 

03039 

338 

7.51 

0.02326 

0.1677E+07 

0.0069 

77.15 

03040 

3.21 

734 

20 


Table  4. 


Initial  state  of  a  Halon  1301  (CBrF3)  jet  from  Figure>l>type  vessel  with  characteris¬ 
tics  of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  rj  =  lOrjj  = 
lOCn^rf,  =  0.07378in. 


Pal  ~  ISOOkg/m^ 
=  <>.17kg/in^ 
V  =  215K 


t 

[S] 

^DV 

[Pa] 

V 

[m/s] 

*1 

[kg/s] 

dMj^/dt 

[kg/s] 

0 

0.4135E+07 

1.0000 

142.83 

0.5181 

9.95 

9.26 

0.00204 

03687E+07 

0.9030 

135.28 

0.5205 

9.43 

8.69 

0.00407 

0J338E+07 

0.8111 

128.99 

0.5224 

8.99 

8.22 

0.00611 

0.3057E+07 

0.7236 

123.64 

0.5239 

8.62 

7.83 

0.00814 

0.2826E+07 

0.6397 

118.99 

0.5252 

8.29 

7.50 

0.01018 

0.2631E+07 

0.5590 

114.91 

0.5263 

8.01 

731 

0.01222 

0.2465E-I-07 

0.4812 

11138 

0.5273 

7.76 

6.95 

0.01425 

0.2322E+07 

0.4058 

108.03 

0.5281 

7.53 

6.73 

0.01629 

0.2197E+07 

03327 

105.08 

0.5289 

732 

6.52 

0.01832 

0.2086E+07 

0.2616 

102.40 

0.5295 

7.14 

634 

0.02036 

0.1988E+07 

0.1924 

99.94 

0.5301 

7.97 

6.17 

0.02240 

0.1899E+07 

0.1248 

97.67 

0.5307 

6.81 

6.02 

0.02443 

0.1820E+07 

0.0588 

95.58 

0.5311 

6.66 

538 

Table  5. 


Initial  state  of  a  Freon  22  (CHCIF2)  jet  from  Figure- 1-type  vessel  with  characteris¬ 
tics  of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  r,  =  20rT  t  = 
20Co^%  =  0.1476m.  ^ 


Pal  =  1320kg/m^ 
Plag  =  3.58kg/m^ 
T^a  =  232K 


t 

[S] 

^DV 

[Pa] 

Vu 

[m/s] 

[kg/s] 

dM^^/dt 

[kg/s] 

0.00000 

0.4135E+07 

1.0000 

32.21 

03172 

536 

11.53 

0.00179 

0J663E-I-07 

0.9030 

30.13 

03158 

5.01 

10.86 

0.00358 

0J299E+07 

0.8117 

28.45 

03148 

4.73 

1030 

0.00537 

0J008E-i-07 

0.7248 

27.05 

03139 

4.50 

9.83 

0.00716 

0.2769E-I-07 

0.6417 

25.86 

03132 

430 

9.43 

0.00895 

0.2570E-I-07 

0.5620 

24.83 

03126 

4.13 

9.08 

0.01073 

OJ401E+07 

0.4851 

23.92 

03121 

3.98 

8.77 

0.01252 

0.2255E-t-07 

0.4108 

23.12 

03117 

3.85 

8.49 

0.01431 

0.2128E-i-07 

03389 

22.40 

03113 

3.73 

8.24 

0.01610 

0.2016E-I-07 

03690 

21.75 

03109 

3.62 

8.02 

0.01789 

0.1917E+07 

0.2010 

21.16 

03107 

3.52 

7.81 

0.01968 

0.1829E+07 

0.1347 

20.62 

03104 

3.43 

7.62 

0.02147 

0.1749E-I-07 

0.0701 

20.12 

03102 

335 

7.44 

0.02326 

0.1677E+07 

0.0069 

19.67 

03099 

3.27 

7.28 

22 


Table  6. 


Initial  state  of  a  Halon  1301  (CBrFj)  jet  from  Figure>l-type  vessel  with  characteris¬ 
tics  of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  r,  =  20r,  t  = 
20CB“rN .  0.1476m.  I  U 


—  ISOOkg/m^ 
Pug  =  6-17kg/m^ 
Tu  =  215K 


t 

[S] 

^DV 

[Pa] 

1 

Vu 

[m/s] 

[kg/s] 

dMj^dt 

[kg/s] 

0.00000 

0.4135E+07 

1.0000 

38.48 

0.5583 

10.73 

8.49 

0.00204 

0J687E+07 

0.9030 

36.17 

0.5566 

10.08 

8.03 

0.00407 

0J338E+07 

0.8111 

34.28 

0.5552 

9.56 

7.66 

0.00611 

0J057E+07 

0.7236 

32.69 

0.5541 

9.11 

733 

0.00814 

0.2826E+07 

0.6397 

3133 

0.5532 

8.73 

7.05 

0.01018 

0.2631E+07 

0.5590 

30.15 

0.5524 

8.41 

6.81 

0.01222 

0.2465E+07 

0.4812 

29.11 

0.5518 

8.12 

6.59 

0.01425 

0.2322E+07 

0.4058 

28.19 

0.5512 

7.86 

6.40 

0.01629 

0.2197E+07 

0J327 

2736 

0.5507 

7.63 

6.22 

0.01832 

0.2086E+07 

0.2616 

26.60 

0.5503 

7.42 

6.06 

0.02036 

0.1988E+07 

0.1924 

25.92 

0.5499 

7.23 

5.91 

0.02240 

0.1899E-I-07 

0.1248 

2539 

0.5496 

7.05 

5.78 

0.02443 

0.1820E+07 

0.0588 

24.71 

0.5492 

6.89 

5.65 
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p 


ftth  in  the  enthalpy  -  pressure  plane  of  the  thermodynamic  state  of  the  initially 
liquid  agent  as  it  flows  from  inside  to  outside  the  pressure  vessel  and  achieves  a 
stable  two-phase  state  at  the  inital  section. 
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AT  z=Zj : 


I, A 


-(dMj  ^e/dt)/(itrj) 


MASS  FLUX  OF  GASEOUS  AGENT 


-(dMj  ^/dt)/()tr*) 


MASS  FLUX  OF  LIQUID  AGENT 


Figure  4.  Sketch  of  the  velocity  and  mass  fluxes  of  agent  gas  and  agent  liquid  droplets  at  the 
initial  section,  z  —  Zj. 
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[Bd]  ''^d 


t[s] 


Figure  5. 


Plot  of  Table  3  results  for  P^y  as  a  function  of  time  during  vessel  discharge  of 
I^%on  22. 


0  0.005  0.01  0.015  0.02  0.025 


Figure  6. 


t[s] 


Plot  of  Table  3  results  for  V|^  as  a  function  of  time  during  vessel  discharge  of 
Freon  22. 


t[s] 


Figure  7. 


Plots  of  IJible  3  results  for  C  and  Xj  as  functions  of  time  during  vessel  discharge  of 
Rreon  22.  ^ 


30 


dM/dt  [kg/s] 
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Figure  8. 


t[s] 


Plots  of  Ihble  3  results  for  dMj^^dt  and  dMj^o/dt  as  functions  of  time  during 
vessel  discharge  of  F^n  22. 
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